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, Abstract. The initial value problem for the conservation law dtu + (— A)"/^m + V • 

' /(^^) = is studied for a £ (1,2) and under natural polynomial growth conditions 

imposed on the nonlinearity. We find the asymptotic expansion as oo of 

solutions to this equation corresponding to initial conditions, decaying sufficiently 
fast at infinity. 



^ , 

p^! 1. Introduction 

I We study properties of solutions to the initial value problem for the multidimensional 

^ ■ conservation law with the anomalous diffusion 
» . 

^ - (1.1) + (-A)"/2y + V ■ /(m) = 0, xeM^t>0, 

(1.2) m(x,0)= Mo- 

Here, we always impose the standing assumption 1 < a < 2. Moreover, we assume 
that the C^- vector field f{u) = (/i(m), . . . , fd{u)) is of a polynomial growth, namely, 
00 ! it satisfies the usual estimates 

^' (1-3) \f{u)\<C\u\'' and \f{u)-f{v)\<C\u-v\{\u\''-^ + \v\''-^) 

for some constants C > 0, g > 1 and for all m, f G M (in fact, assumption (11. 3p can be 
slightly relaxed in some parts of our considerations, cf. Remark 12. 3[ below). 
O ■ Linear evolution problems involving fractional Laplacian describing the anomalous 

diffusion (or a-stable Levy diffusion) have been extensively studied in the mathemat- 
ical and physical literature (see, e.g., [H]). The probabilistic interpretation of non- 
linear evolution problems with an anomalous diffusion, obtained recently by Jourdain, 
Meleard, and Woyczyhski [9], motivated us to study fll.ll) - fll.2p . The authors of [9] con- 
sidered a class of nonlinear integro-differential equations involving a fractional power of 
the Laplacian and a nonlocal quadratic nonlinearity represented by a singular integral 
operator. They associated with the equation a nonlinear singular diffusion and proved 
propagation of chaos to the law of this diffusion for the related interacting particle 
systems. In particular, due to the probabihstic origin of fll.ll) - fll.2p . the function u{-,t) 
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should be interpreted as the density of a probability distribution for every t > 0, if the 
initial datum is so. 

Although, the motivation for this paper comes from the probability theory, our 
methods are purely analytic. Hence, if X{t) is the symmetric a-stable Levy process, 
its density of the probability distribution pa{x,t) is the fundamental solution of the 
linear equation 

(1.4) dtv + (-A)"/2^; = 0, 

hence, Pa can be computed via the Fourier transform Pa{^,t) = e"*'^'". In particular, 

Pa{x,t) = r°'/"P„(a;r^/"), 

where Pa is the inverse Fourier transform of e~'^'° (see [8], Ch. 3] for more details). It 
is well known that for every a G (0, 2) the function is smooth, nonnegative, and 
satisfies the estimates 

(1.5) < Paix) < C{1 + |x|)-(°+'^) and |VP«(x)| < C(l + |x|)-("+'^+^) 
for a constant C and all x G M*^. Moreover, 

(1.6) Pa{x) = co|a;|-("+'^) + O , as |x| ^ oo, 
and 

(1.7) VP„(X) = -Ci ^ Q (^|^|-(2a+d+l)^ ^ ^ 

where 

CO = a2-^vr-('^+^)/^in(W2)r(^)r(|), 

and 

ci = 27ra2"-V-('^+^)/2sin(a7r/2)r(^±|i^)r(|). 

We refer to [3] for a proof of the formula (11. 6p with the explicit constant Cq. The 
optimality of the estimate of the lower order term in (11. 6p is due Kolokoltsov [TUl Eq. 
(2.13)], where higher order expansions of Pa are also computed. The proof of the 
asymptotic expression (11.71) and the value of Ci can be deduced from (II. 6p using an 
identity by Bogdan and Jakubowski [U Eq. (11)]. 

The asymptotic formula (II. 6p for the kernel Pa plays an important role in the theory 
of a-stable processes. The main goal in this work is to present a method which allows to 
derive analogous asymptotic expansions as |a;| ^ cxd of solutions to the Cauchy problem 
(ll.ip - (ll.2l) . In the next section, we recall several properties of solutions to (II. 10 - 01. 2p 
and we state our main results: Theorems 12.11 and 12.41 In Section |3l we gather technical 
space-time estimates of solutions to (ll.ip - (ll.2l) . The proofs of Theorems 12.11 and 12.41 
are contained in Section HI 

Notation. The L^-norm of a Lebesgue measurable, real-valued function v defined 
on M'' is denoted by H^Hp. In the following, we use the weighted L°° space 

(1.8) = {v e L^^iW^) : \\v\\L^=esssnp^^^,\v{x)\{l + \x\f<oo}, 
for any > 0, and its homogeneous counterpart 

= {t; G (R^{0}) : ^ esssup,,j,.Kx)||xr < oo}. 
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The constants (always independent of x) will be denoted by the same letter C, even if 
they may vary from line to line. Sometimes, we write, e.g., C = C{T) when we want 
to emphasize the dependence of C on a parameter T. 

2. Main results 

It is well known (see [H [5l [6]) that given uq G L^(R'^) and 1 < a < 2, the initial 
value problem fll.ll) -( fL2l) has the unique solution u G C{[0,oo),L\R'^)). Moreover, 
this solution satisfies 

u G C{{0,oo),W^'P{R'^)) 
for every p G [1, oo] and the following estimates hold true (see \i, Cor. 3.2]) 

(2.1) Mt)\\p<cr^^'--^^\\uo\\i 

for all t > and C independent of t and of Uq. Under the additional assumption 
Mo G LP(M.'^), the corresponding solution satisfies u G C([0, oo), Lp(R°')) together with 
the estimate 



(2.2) \Ht)\\p<\\uo\ 



V 



Below, in Proposition 13.31 we complete these preliminary results providing the es- 
timates of solutions to fll.ll) - fll.2p in weighted L°°-spaces. In particular, if uq G L'^j^^ 
(cf. dLl])), then the corresponding solution of ([III])-([L2]) satisfies u G C([0,T],L^+J 
for every T > 0. Such a result was already obtained in the one dimensional case, see 
[3, Sect. 2]. In Section [3l we state and prove its multidimensional counterpart for 
the completeness of the exposition. We complement this result with additional esti- 
mates for the gradient of the solution, which will be useful in the proofs of asymptotic 
formulas in Section HI 

Let us recall that, when studying the large time behavior of solutions for the prob- 
lem f ll.ip - fll.2l) . an important role is played by the critical exponent 

^ a — \ 

Indeed, using the terminology of [2j the behavior of solutions as t ^ oo is genuinely 
non-linear when g = g, is weakly non-linear when q > q, and is (expected to be) 
hyperbolic when 1 < q <q. 

In this paper, in the supercritical case g > g, as well as for g = g provided ||mo||i 
is sufficiently small, we will improve the space-time estimates of [TJ Sect. 2], showing 
that 

(2.3) \uix,t)\<Cpaix,l+t), 

for all X G M'^, t > 0, and C > independent of x,t. Under the additional assumption 
that Vmq G L'^^^_^_l, we will also prove that 

(2.4) \\Vuit)\\L^^^^^<Cil+t), 



see Theorem 13. 5[ below. In other words, Vm(x, t) has the same space-time decay profiles 
as Vpa{x, I + t) (cf. the second inequality in (11.51) ). 

Furthermore, we make evidence of the second critical exponent, namely, 

a + a 
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playing an important role in the study of the pointwise behavior of solutions as |a;| — > 
oo. The following theorem explains the role of g*, showing that any decaying solution 
has a precise spatial asymptotic profile. Here, we denote by Sa(t)uo{x) = Paif) *uq{x) 
the solution of the linear equation (11.41) supplemented with the initial datum uq. 

Theorem 2.1. Assume that a G (1, 2), and let u = u{x, t) be the solution of (II. II) -( 11121) 

with the nonlinearity satisfying (li.5| . and with Uq G L"^^^. 
(z) Then, for allt>0, x e M.'^, 



u{x, t) = Sait)uo{x)+j^^^^^^ ■ j^j (t - s)f{u{y, s)) dyds 



as \x\ — > oo, 



(2.5) 

+0 (^max{ |x| ; |x| ^("+"^+2) }) , 

uniformly in any time interval t G [0,T], T > 0. This conclusion is interesting only 
when the last term on the right hand side of (1^.51) is the lower order term as \x\ —>■ oo: 
this happens when q > q* . 

(a) The conclusion at the point (i) can he improved under the additional assumption 
Vmo G -L^rf+i, replacing the third term on the right hand side of ( \2.5\^ by 



as \x\ — > oo. 



Now, this conclusion is interesting also for 1 < q < q* . 

(Hi) If u satisfies inequality l \2.3\} for all x G M."^, t > 0, and C > independent 
of X, t, then one can make precise the behavior for large t of the remainder term in 
relation {\2.5\i . replacing it by 



0((1 + t)^max{|x|-^(°+'^) ; }) , 



as \x\ oo, 



uniformly in t & [0, oo), for some exponent N = N{a, q, d) < 3, independent on uq. 

If, in addition, the solution satisfies inequality i \2.4^ , the conclusion at the point (ii) 
can be improved replacing the remainder term by 

0((1 + t)^max{|xr'^("+^)-i ; }) 

for some exponent N = N{a, q, d) < 3, independent on uq, and the convergence as 
\x\ — > oo holds true uniformly in t E [0, oo). 

It follows from the Duhamel formula that the solution of the Cauchy problem (II. ip - 
(II. 2p satisfies the integral equation 

(2.6) u{t) = S^{t)uo - [\s^{t-s)-f{u){s)ds. 



Hence, it is possible to give a heuristic explanation of the role of space-critical exponent 
q = q*, simply, by looking at the integrand of the second term on the right hand side 
of (12. 6p . Indeed, the kernel of VS'^ behaves as as |x| oo (cf. (I1.6P ). 

whereas \ f{u{x,t))\ < C\x\~''^'^~^°'^ for u{t) G L'^^^- Then, it is natural to expect that 
the large space asymptotics is influenced by the competition between these two decay 
rates as |x| — > oo. In fact, the proof of Theorem 12.11 (given in Section H]) consists in 
finding the asymptotic expansion of the second term on the right hand side of (12. 6p 
and the equality between these two decay rates occurs precisely when q = q*. 
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Remark 2.2. It is worth observing that this type of asymptotic expansion of solutions to 
convection-diffusion equations is specific of the fractional nature of the diffusion opera- 
tor (— A)"/^ and is caused by the algebraic decay of the fundamental solution pa{x,t). 
For the viscous Burgers equation, or for multidimensional diffusion-convection equa- 
tions with standard dissipation {i.e., with the usual Laplacian) Theorem 12.11 remains 
valid, but it is not interesting because the coefficient ci vanishes in the limit case a = 2. 

Remark 2.3. The conclusion (i) of Theorem 12.11 remains valid under more general 
assumptions on the nonlinearity. What we really need is that / is a C^-vector field 
such that < c{R)\u\'^ for some g > 1, a continuous nondecreasing function c(-) 

on [0, oo), and all \u\ < R. For the part (ii), we need also a similar condition for 
/', namely, \f'{u)\ < Ci{R)\u\'^~^ for \u\ < R. On the other hand, the present form 
of Theorem I2.1[ iii is no longer valid for such more general nonlinearities. Our more 
stringent assumption fll.3p allows us to present the essential ideas avoiding uninteresting 
technicalities in the proofs, in particular, separating the cases of large and small u in 
our estimates. Moreover, such an assumption is well suited for studying self-similar 
solutions. 

For the homogeneous nonlinear term V ■ f{u) = b ■ V(n|n|'^^-'^) with a fixed 6 G M"' 
and with the time-critical exponent q = q, the authors of [2j constructed a family 
of self-similar solutions um = UM{x,t) of equation (11.11) . Those functions satisfy the 
scaling relation 

(2.7) MM(a;,t) =t'^/"t/A/(xri/") where Um{x) = um{x,1) 

for all X G M*^ and t > 0. Moreover, each of them is the unique solution of the initial 
value problem 

(2.8) dtu + {-A)"/\ + b ■ V(m|m|("-i)/'^) = 

(2.9) u{x, 0) = MSo 

for a G (1,2) and M > 0, where 6o is the Dirac delta. We refer the reader to [2] for 
more information concerning solutions of problem (I2.8p - (I2.9I) . 

In this paper, we complete results from [2] providing space-time estimates of those 
self-similar solutions. First, in Corollary 13.61 below, we establish, for sufficiently small 
M > 0, the estimate 

(2.10) < um{x, t) < Cp^{x, t) for all x G M'^ and t > 0, 

We conjecture that such estimate remains true without the smallness assumption im- 
posed on M. Inequality fl2.10p plays a crucial role in the proof of the following asymp- 
totic expansion of the self-similar kernel Um- 

Theorem 2.4. Assume that 1 < a < 2 and q> q* . Let um be a self-similar solution 
of (l^.gp - O^.gp . satisfying the estimate (1^. iOp . Then the self-similar profile Um{,x) = 
um{,x, 1) has the following behavior as \x\ oo: 

(2.11) f/^,(x) = MP„(x) + ^||f/M|||^p^ + 0(max{|x|*^^ 

The asymptotic expansion of solutions to (II. ip stated in (12. 5p and in (12. lip can be 
viewed as the true counterparts of the well-known result for the a-stable distribution 
recalled in (11.61) . 
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3. Preliminary space-time estimates 
We begin this section by the study of the solution of the hnear problem 

(3.1) dtv + {-A)''^\ = 0, v{x,0) = vo 
denoted by 

V{X, t) = Sait)Vo{x) = Pai-, t) * Vo{x). 

The following lemma contains a direct generalization to of estimates from [TJ Lemma 
1.40]. By this reason, we sketch its proof only. 

Lemma 3.1. Assume that Vq G L'^^^- There exists C > independent of Vq and t 
such that 

(3.2) \\S^{t)vo\\oo < Cmin{r^/"||t;o||i,||t^o||oo}, 

(3.3) \\Sa{t)vo\\L^^^ < Cil + t)\\vo\\L^^^, 

(3.4) l|V5Jt)t;o||L-, < Ct-'/''\\vo\\L^^^ + Ct'-'/''\\vo\\i, 

Proof. Estimate (13. 2p results immediately from the Young inequality applied to the 
convolution Sa(t)vo = Paif) * Vq, due to the identities 

|K(t)||i = l, |K(t)||oo = r'^/l|Pa||oo foralH>0. 

Since |fo(a;)| < C{l + \x\)^^'^^"\ by the asymptotic properties of the kernel Pq,(x, 1) = 
Pa{x) (cf. (11.61) ). we immediately obtain |fo(a;)| < Cpa{x,l) for all x E M.^ and 
a constant C > independent of x. Consequently, by the semigroup property, we 
conclude 

\\Sait)vo\\L^^, < C\\S^it)p^il)U^^^ = ciM + 1)|U^^^ < C(l + t). 

Now, replacing vq by t'o/||fo||L^^^ we obtain (13.31) . 
To prove (13. 4p . we use the pointwise estimate 

(1 + ixD^^'^ < c(i + ii/i)"+^ + c\x - yr+^ 

valid for all x,y eM. and a constant C > 0, and we apply the Young inequality. We 
get 

\\VS^it)vo\\L^^^ < C||Vp«(t)||i||.;olU^^, + C||Vp„(t)||^^^J|i;o||i 

and (13.41) immediately follows. □ 

Under an additional information on the gradient of vq, we can obtain analogous 
estimates for VSait)vo. In order to give a precise statement, let us introduce the space 

(3.5) E^^, ^{ve <f (M'^) : ||.;||^^^, = \\vh^^^ + || Vt;|U^^^^^^ < oo}. 

Lemma 3.2. Assume that vq G Ea+d- There exists C > independent of Vq and t 
such that 

(3.6) ||V5,(t)t;o||oo < Cmin{r('^+i)/°||t;o||i;ri/"||t;o||oo; ||Vt;o||oo}, 

(3.7) \\Sa{t)vo\\E^^, < C{l + t)\\vo\\E^^,, 

(3.8) ||V5,(t)i;o|U.+, < Ct-'^^voWE^^. + Ct'-'/^voh 
for all t > 0. 
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Proof. Estimate (13.61) is the straightforward apphcation of the L^-L°° convolution in- 
equahties. In order to prove fl3.7p using the radial symmetry of Pa{-,t), we see that, for 
all R> 0, f^^ ^p{y, t) dy = 0, where Br denotes the ball centered at the origin and of 
radius R. Hence, 



VSait)vo{x)= [voix-y)-Vo{x)]Vpaiy,t)dy+ vq{x -y)Vpo,{y,t) dy. 

-'l?/l<l^l/2 J\y\>\An 

This decomposition shows that, for some constant C > 0, the quantity \V So,{t)vQ{x)\ 
can be bounded from above by 

C|x|-("+'^+i)||V^;o||l- / |i/||Vp(2/,t)Mi/ + a|x|-(°+'^+i) / W{y)\dy, 

which implies 

(3.9) \\VSo.it)voh^ <c(\\Vvoh^ +tp^^^^ 

Now, estimate (13. 7p follows from (13. 3p . (13. 6p and from the bound for the homogeneous 
norm (13. 9p . 

Let us prove (13. 8p . By (13. 4p and the inequality 

\\V^S^{t)vo\\oo < ||V5,(t)||i||Vt;o||oo < Ct-i/"||Vt;o||oo, 
we see that we only have to establish the following estimate in the homogeneous space 

r oo 

(3.10) \\V''S^{t)vo\\L^ <Ct-'^''\\vo\\E^^, + Cf-'/''\\voh. 
To prove (13.100 . we consider the decomposition 

V^Sa{t)vo{x) = (Ji + J2 + J3){x,t), 

where 



Ji(x,t)= / [vo{x -y) -vo{x)]V pa{y,t)dy, 

'l2/l<k'l/2 



J2(x,t)= / vo{x -y)V Pa{y,t)dy, 

J3{x,t) = -Vo{x) V^pa{y,t)dy 

(note that /jg^ "V^Paiyyt) dy = 0). /^From the well known estimate (see [TO] ) 
(3.11) \V^Pa{x)\ < C{1 + |x|)-("+'^+'), 

we deduce J^a \y\ I V^Pq,(?/, t)\dy < Ct~^/°'. Then, the apphcation of the Taylor formula 
in the integral defining Ji yields 

|Ji(x,t)| < Cri/"|x|-(°+'^+i)||Vi;o|U- 

To deal with the terms J2 and J3, we use two different pointwise estimates of 
VVq,(x, t) resulting from (13. lip : 
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and 



which imply 



|^2(a;,t)| < sup \V^Pa{y,t)\ \vo{x-y)\dy 

\y\>\x\/2 J\y\>\x\/2 



and 



\J3{x,t)\ < C|a;r(-+^)||^o|U2=+. / \V'p4y,t)\dy 

-'|j/|>kl/2 

II II U||1J^_^^ 

Combining all these inequalities yields (I3.10p . □ 

We are in a position to construct solutions of the Cauchy problem ( ll.ip - fll.2l) in the 
weighted space L"^^^. 

Proposition 3.3. (i) Let a G (1,2) and q > 1. Assume that u is a solution of the 
Cauchy problem fll.ll) - fll.2l) with the nonlinearity satisfying (11 ■3p . If Uq G L"^^^, then 

(3.12) u G C{[0, T], L^+rf) for each T > 0. 

(a) Under the more stringent assumption uq G Ea+d, cf. (1 5*. 51) . we have also 

(3.13) u G L°°([0, T], Ea+d) for each T > 0. 

Proof. In order to prove (13.121) . it suffices to show that the nonlinear operator 

T{u){t) = Sa{t)uo - [ VSa{t-T)f{u{T))dT 

Jo 

has the fixed point in the space 

Xr = {wGC([0,T],L-,) : sup \\uit)U^ < 

te[o,T] 

As usual, we work in the ball B{0,R) = {n G C([0,T],L^^) : sup^gjQ j.] ||^^(i) < 
R}, where R = MUmoIU^,^^ and M > is a large constant, and T > 0. Combining 
inequahty (13. 4p with assumption (11.31) we get 

\\vSa{t)fiu)U^^^ < cr^/llkniL^^^ + cti-^/lkll^ 

(3.14) < Ct-y-il + t)\\ur^'\\uh^^^. 

Applying now inequalities (l3.3p -( l3rT4l) we can estimate, for u G B{0,R), 

\\Tium\L^^^ < C{l + t)\\uo\\L^^^ 

-t 



+CR^-' [ {t-rr'/-{l + {t-T))\\u{r)U^^^dT 
Jo 

< R/2 + CM''-^\\uo\\tJ Rt^'^'^'il + t) 



a + d 

< R, 



provided that < t < T and 



T<Cmin{l,||Mo||i°^'^"'^/^"^'^}, 
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with C > small enough. 

In the same way, for all u,u & B{0, R), 

Wnum - Tm)h^^^, < CW-' [\t - r)-V"(l + (t _ rmuir) - uir)\\,^^^ dr. 



Jo 

The Banach fixed point theorem now guarantees the existence of a local-in-time solu- 
tion. In the next step, such solution must be extended globally-in-time. The argument 
is standard: we fix T > arbitrarily large and using that ||M(t)||oo < C on [0,T] (see 
inequality (12.21) ). we show that II "^^(1^) II L^+d ^'^^^ blow up on [0,T]. Indeed for some 
constants Ci, C2, ■ ■ ■ , depending on T, for < t < T we have 

\\u{t)h- < Ci + [\t - rr'/-\Hr)U^^^ dr. 



Jo 

Iterating this inequality and applying Fubini's theorem we get 



<c, + c, [\t~Ty-'/^u{T)u^^^dT. 

Jo 



We repeat this argument until we obtain the integrand factor [t — r) with a positive 
exponent; here, only a finite number if iterations are needed, since a > 1. This leads 
to \\u{t)\\L^^^ < C5 + Ce /o \\uir)\\L^^^dT and finally to \\u{t)\\L^^^ < C^expiCet) by 
the classical Gronwall lemma. 

To prove of (I3.13P under the stronger assumption uq G Ea+d, one could proceed in the 
same way, replacing the space L'^^^ with Ea+d (and using the estimates of Lemma [321) • 
However, this argument would require additional restrictions, such as inequalities of 
the form \f'{u) — /'(f) | < C\u — t'|(|'u|''~^ + |f |''~^), which are not fulfilled for some 
nonlinearities satisfying (II. 3p with q < 2. 

Let us proceed in a slightly different way. First of all we have, by [2l [6], Vu{t) G 
L°°([0, T], L°°{R'^)) for all T > 0. We rewrite the integral equation ([M]) in the following 
way 

Vu{x,t) =VSa{t)uo{x) 

(3.15) 



+ / \Vpaix-y,t- s)Vf{u{y,s))dyds. 

y\<\^\f2 J\y\>W2y 

It follows from condition (11.31) that \f'{u)\ < C|n|'^^^, hence, for every u satisfying 
fl3:T2|) we have 

(3.16) \VfHy,s))\ < C(l + |y|)-(''-i)("+'^)|V«(y,.)| < Cil + |^|)-(^~i)("+^), 

for a positive constant C = C{T) and all y G M"', s G [0,T]. Combining (13.161) with 
(13. 7p and with the decay estimate \Vpa{x,t)\ < Ct\x\^^"'^'^^^\ we get from (I3.15p the 



preliminary inequality 

(3.17) \\/u{x,t)\ < C(l + |x|)-("+'^+^) +C(1 + |x|)-(°+'^+^)+'^^ +C(1 + |x|)~('^-^)(°+'^) 

for some constant C = C(T) > 0, all x G M'^, t G [0, T], and with qi = d. Since q > 1, 
now we can use this inequality to improve the estimate in (I3.16p . This allows us to 
replace qi with some < q2 < qi and to improve also the estimate of the third term 
in (I3.17p . After finitely many iterations of this argument (more and more iterations are 
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needed when q approaches 1), we get \Vu{x,t)\ < C(T)(1 + for all x 

andtG[0,T]. □ 

Let us now recall a singular version of the Gronwall lemma. This fact seems to 
be well-known, we state it, however, in the form which is the most suitable for our 
application and we prove it for the completeness of the exposition. 

Lemma 3.4. Assume that a nonnegative and locally hounded function h = h{t) satisfies 
the inequality 

(3.18) h{t) <Ci{l + t) + C2 [ {t- r)-'^(l + r)-'h{T) dr 

Jo 

for some a G (0, 1), 6 > 0, positive constants Ci and C2, and all t > 0. If a + b > 1, 
then h{t) < C(l + 1) for all t > and C independent of t. The same conclusion holds 
true in the limit case a + h = 1 under the weaker assumption 

(3.19) h{t) <Ci{l + t) + C2 [ {t- r)-V-^/i(r) dr 

Jo 

provided C2 is sufficiently small. 

Proof. If a + 6 = 1, we deduce from fl3.19p the following inequality 

h{t) <Ci{l + t) + C2K{a, b) sup /i(r), 

0<T<t 

where 

K{a,h)= I {t-T)-^T-^ dr = I {I - s)-^s-^ ds. 
Jo Jo 

Consequently, supo<T-<j 

^(^) ^ i-cfx(a,b) (l + ^) provided C2 < l/K{a, h). 
In the case a + h> 1, using (I3.18p . we write b = bi+rj with a + bi = 1 and r/ > 0, and 
we fix ti > such that 6*2(1 + tiy < 1/K{a, bi). Now, splitting the integral in (I3.18p 
at ti yields 

h{t) < C{1 +t) + C2K{a, 6i)(l + ti)-'' sup h{T) 

0<T<t 

for some C > independent of t. The conclusion of Lemma [3.41 now follows. □ 

If the exponent q in the assumptions on the nonlinearity f 1 1.3 1) is larger than the 
time-critical value q, we can improve the space decay estimates from Proposition 13.31 
through the following space-time decay result. 

Theorem 3.5. (i) Let a G (1,2). Assume that u = u{x,t) is a solution of the Cauchy 
problem fli.il) - fli.^) . where the nonlinearity f satisfies (\1.3h with q>q = l + {a — l)/d 
and Uq G L'^^^. There exists C > (depending on uq hut independent ofx,t) such that 

(3.20) \u{x, t) I < Cpa{x, 1 + t) for all x G M"' and t > 0. 

The same conclusion holds true for q = q provided ||mo||i is sufficiently small, 
(a) Under the more stringent assumption Uq G E^+d we have also 

(3.21) ||Vn(t)|U^^^^ <C(l+t). 
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Proof. First recall that by estimates (12. ip and fl2.2p with p = oo, the solution satisfies 

||«(t)||oo<C(l + t)-'^/". 

Hence, to establish (13.200 . it suffices to prove 

(3.22) \\uit)U^^^<Cil + t). 
Indeed, the inequahty 

^(x,t)^min|(l + t)-^/^ ^-^±1-^1 <Cp,(x,t + l). 

is the consequence of the elementary estimate 

gix,t) < (l + t)-'^/"min{l; |x(l + t)"^/"!""-'^} 

and the asymptotic formula (11.61) (implying, in particular, that min{l ; |a:|~"^'^} < 
CPaix) for all X e M'^ and a constant C > 0). 

In the proof of (I3.22p . we use the integral equation (12. 6p . hence we begin by the 
preliminary estimate (resulting from (13.40 and from the hypothesis (11.30 ) 

||V5„(t-r)/(n(r))|U^^^ < 0(1 - r)-'/-\\u{r)\\l^'\\u{T)U^^^ 

+Cit-ry-'/-\\uir)\\l. 

Moreover, since by (12.11) and (12. 2p with p = q, the solution satisfies the decay estimate 

(3.23) ||^i(r)||^<C(l + r)-^(''-i)/", 
we have the following inequalities 

{t-Ty-^/''\\u{TWdT<C [ (t-r)i-i/"(l + r)-'^(''-i)/" dT<C{l + t) 



which are valid for 1/a + d{q — l)/a; > 1. 

Consequently, after computing the L^^-norm of equation (12. 6p and using estimate 
(13. 3p we arrive at 

(3.24) \\umL- <C{l + t)+C [\t - r)-i/"(l + r)-'^(^-i)/"||n(r)|Uj» ^ dr. 



Jo 

In the time-critical case l/a + d{q — l)/a = 1 (i.e. for g = g) we proceed analogously, 
however, now we use the estimate 

(3.25) h(r)||oo<Cr-^/ll«o||i 

with a constant C independent of uq and t. Hence, we obtain the following counterpart 
of inequality (I3.24p 

(3.26) \Ht)h^^^ < C{1 + t) + CWuoWr' fit - r)-i/V-^(«-i)/"||n(r)|U^.^^ dr. 



Finally, the singular Gronwall lemma (Lemma 13. 4p applied to inequalities (I3.24p and 
(I3:26|) completes the proof of ^TM) . 

To prove inequality (13.211) one should follow exactly the same argument as for the 
proof of (I3.22p . putting everywhere E'^+d-norms instead of the corresponding i^^+d" 
norms, and applying Lemma 13.21 instead of Lemma 13. 1[ □ 
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We conclude this section with estimates of self-similar solutions to problem fl2.8p - 

Corollary 3.6. // the constant M > in l\2.!M is sufficiently small, then the corre- 
sponding solution of problem i\2.8^ - (\2.9^ satisfies 

(3.27) < um{x, t) < Cpa{x, t), for all xeR'^,t>0, 
with C = C{M, a,d) > independent of x and t. 

Proof. Let us recall that the solution of fl2.8l) - fl2.9p has been constructed in p] as the 
limit of the rescaled functions u^{x,t) = X'^u{Xx, X^'t), where u = u{x,t) is the fixed 
solution of equation (12.81) supplemented with the nonnegative initial datum m(-,0) = 
uq e C^(M'^) such that J^^uoix) dx = M. By Theorem 13. 5[ used in the critical case 
q = q, the rescaled family u'^ satisfies 

(3.28) \u^{x, t) I < CX'^PaiXx, 1 + A"t) = Cp„(x, A"" + t) 

for all X e W^, t > 0, and a constant C = C{M,a,d) independent of x,t,X, provided 
M > is sufficiently small. Since u^{x,t) UM{x,t) as A — oo almost every where 
in {x,t) (see p] Lemma 3.7]), passing to the limit in (13.281) we complete the proof of 
estimate (13:271) . □ 



4. Asymptotic profiles 

In this section, we derive the asymptotic expansions from Theorems 12.11 and 12 . 4[ Let 
us recall that all positive constants, which appear here, are independent of x and t and 
are denoted by the same letter C. 

Proof of Theorem \2.1\ Let us consider the nonlinear term appearing in the integral 
equation (12. 6p . 

J^{u){t) = j^j VMdPa(a; -y,t- s)f{u{y, s)){s) ds. 

In order to find an asymptotics of A/" for large we define two remainder functions 
7l{x,t) and 7li{x,t), through the relations 



ft 

N'{u){x,t) = I I f{u{y,s))Vpa{x,t - s)dyds + ni{x,t) 



(4.1) 



Ci X 



\x 



a+d+2 



t 

{t — s)f{u{y, s)) dyds — lZ{x, t). 



Here, Ci is the constant appearing in relation (II. 7p . Hence, it follows from the integral 
equation (12. 6p that 

(4.2) u{x,t) = S^{t)uo{x) + -^^^ f [ {t-s)f{u{y,s))dyds + n{x,t) 
and it remains to estimate TZ{x,t). 
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Computing the difference of the two expressions of N in (14. 11) we deduce a bound 
for 7^ + 7^1, implying 

|7^(a;,t)| < |7^l(x,^)| 



+ C\x\ 



-(2a+d+l) 



JM.'^ 



\f{u{y,s))\dyds. 



Now, the asymptotic expansion (ll.7p . the assumption fll.Sp . and the L'^-estimates f l3.23p 
lead to 

|7^(x,^)|<|7^l(x,t)| + c|xr(2"+^+l) /7 (t-sf 
(4.3) -^o^f 

< |7^l(a;,^)| +C|x|-(2-+rf+i)t2 /" (i + rf^. 



In order to estimate TZi, we decompose it as 7?.i = Ji + ■ • • + I4, where 

Ii{x,t) = / [Vpa{x-y,t- s) -Vpa{x,t- s)]- f{u{y,s))dyds, 

Jo J\y\<\x\/2 

l2{x,t) = - (/ f{u{y,s))dy]Vpa{x,t - s)ds, 

Jo \J\y\>\x\/2 J 

h{x,t) = / / Vpa{x -y,t- s)f{u{y,s))dyds, 

Jo J\y\>\x\/2, \x-y\>\x\/2 

Vpa(2/, t - s)f{u{x - y, s)) dyds. 

l<l^1/2 

In our next two estimates, we use the inequality (which is a consequence of the 
L°°-bound of the solution, see (12.21) ) 

(4.4) \u{y, .)r < C(l + |y|)-("+^ni + 

This leads to 

|/i(x,t)|<qxr("+^+2) f\t-s)[ \y\\u{y,s)\'^dyds 



h{x,t) = [ [ 

Jo J\y\: 



< C|xr("+^+2)t / (l + s)-(''-^)'^/"||M(s)|Uc»^ds. 

Jo "^"^ 

Here, we have applied also the Taylor formula and the bound (13. lip . 
The next two integrals can be bounded by the same quantity, indeed 

|/2(x,t)| + |/3(x,t)| <C|x|-("+^+i) [\t-s) [ \u{y,s)\'^dyds 
(^4g^ -'0 J\y\>\x\/2 

Jo "'^'^ 

The estimate for the last term is 
(4.7) |J4(x,t)| <C|x|-^("+'') r(t-s)-i/"||M(s)||^oo ds. 
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Since we are assuming a > 1, when we compare the exponents of in inequahties (14 .Sp 
and dSD-dMD, we see that 

(4.8) \n{x,t)\ < C max{|x|-'?(°+'^) ; for all |x| > 1 and t G (0,r], 

where C = C (T) > is uniformly bounded with respect to x and t, in any time interval 
t G [0,T]. Part (i) of Theorem 12. II now follows. 

To establish Part (ii), we have only to improve the estimate of the integral (14. 7p . We 
can do it using, in a slightly deeper way, the properties of the fundamental solution 
Pa{x,t). In particular, its radial symmetry implies that 



/ 



'^Pa{y,t- s) ds = 0, 

'l2/l<l^l/2 

SO that 

(4.9) h{x,t)= j / Vpa{,y,t- s)-[f{u{x-y,s))- f{u{x,s))\dyds. 

Jo J\y\<\x\/2 



Owing to the more stringent assumption uq G Ea+d and by Proposition 13. 3[ we deduce 
from the mean value theorem applied to f{u) (recall that \f'{u)\ < C\u\'^~^) 



(4.10) \h\ < c\xr^-^'^'' \Hs)r,ipu{s)u^^^^^ds. 

Replacing inequality (14. 7p with this new estimate shows that the bound (14.81) of the 
remainder term can be improved into 

(4.11) \n{x,t)\ < C max{|xr''("+'^)-^; Ixr^^+'^+^^j for all |x| > 1 and t G (0,T]. 

Hence, Part (ii) of Theorem 12.11 follows. 

Let us prove assertion (iii). When the solution satisfies the additional estimate (12. Sp 
(recall that, by Theorem 13.51 such an estimate holds true at least when either q > q 
or q = g and ||mo||i is small enough), we have ||'w(t)||Lj=^^ < C(l + 1)- In this case, it is 
easy to construct an exponent N = N{a, d, q) such that 

(4.12) |7^(x,^)| < C(l + t)^max{|a;r'?(°+'^); for all |a;| > 1, t > 0. 

Let us explain why < 3. It follows directly from (14. 3 p and from (14. 50 - (14. 7p that 
< max{3 ; g + 1 — 1/a}. However, if g > 2 + 1/a, then we can replace estimate (14. 7p 
with 



(4.13) 



\h{x,t)\ < /\t-s)-l/"||M(s)||^l-2/(a+d)||^^^)||l+2/(a+d)^^ 

Jo 

<C(l + t)3|x|-("+'^+2)_ 

If, moreover, Vm satisfies the additional pointwise estimate ( 12.40 then we can precise 
in a similar way the bound ( 14.110 . Namely, we can replace C = C{T) in ( 14. lip with 
C(l + t)'^. Next, the proof of this claim relies either on inequality (14.100 if 1 < g < 2 
or on the following new estimate of J4 when g > 2 

\h{x,t)\ < C\x\-^-^'^'^ f ||^s)||-i-V("+^)||n(.)||i/i":'^)||VMs)|U. ds. 



The estimates of the other terms remain unchanged. The proof of Theorem (12. ip is 
now complete. □ 
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Proof of Theorem \2.4\ Let mm be a self-similar solution of fl2.8l) - fl2.9p . satisfying esti- 
mate fl2.10p . We consider the integrals /i, I2, I3 and J4 and also the remainder term 
TZ, defined as in the proof of Part (i) of Theorem I2.1[ We treat all these terms pro- 
ceeding as before, but replacing everywhere estimate (14.41) with the estimate (deduced 
from fICTD ) 

(4.14) \uM{y,s)\'^ < Cs-^'^/^P^iy/s'/"), 

with q = q and C = C{M), then making the change of variables y 1— >■ ys^^"" in all the 
space integrals. After some simple computations, we arrive at 



\n{x,t)\ < ct-'^'° 



[2a+d+l) / \x\ \-{a+d+2) / \x\ \ -q{a+d) 
+ ^] + 



Recalling that f{u) = bu^, applying (14.21) to um we get 

UMix,t) = Mp„(x,t) +t^+i/" ■ ^(/ UM{yydy^j^;^^+n{x,t). 
Now, passing to self-similar variables, we deduce that, for all x G M^, 
Um{x) = MP^ix) + ^WUmWIj—^ + 7^M(x), 



where 

UmIx) = 0(^max{|x|-^("+^) ; , as \x\ 00. 

Theorem (12.41) is now established. □ 



Remark 4.1. We conclude observing that the above expression of the remainder term 
TZm{x) can be simplified distinguishing the two cases d = \ and d > 2. Indeed, an 
elementary calculation shows that 

1. In the one dimensional case d = 1 (hence, q = a, and the assumption q > q* 
reads a > a/2), we have 

^ . . / 0(|xr-(-+i)) if v^<a<V3, 11 
^^^(^) = 1 0(N-("+3))' if V3<«<2, 

2. For > 2, it follows 

7^M(x) = 0(|xr^("+'^)) as |x| 00. 



Remark 4.2. Analogously, as in Theorem 12.11 one could remove the restriction q > q* 
from Theorem 12. 4[ provided we have the additional weighted estimate 

(4.15) \\VuMmL^^,^,<Ct. 

We expect that inequality (14.151) can be proved using the scaling argument from the 
proof of Corollary 13.61 below. This reasoning would require, however, some improve- 
ments of estimates from [2]. We skip other details because the goal of this work was to 
present a method of deriving asymptotic expansions of solutions rather than to study 
the most general case. 
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